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We present a theory of the Forster energy transfer between the arrays of donor and acceptor 
molecules lying on the planar metallic mirror. We reveal strong modification of the effective transfer 
rate by the mirror in the incoherent pumping regime. The rate can be either suppressed or enhanced 
depending on the relative positions between acceptor and donor arrays. The strong modification of 
the transfer rate is a collective effect, mediated by the light-induced coupling between the donors; 
it is absent in the single donor model. 
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I. INTRODUCTION 

Forster energy transfer processes are now actively 
studied in various fields that bridge physics, biology and 
medicine. The energy is transferred from the excited 
(donor) system to the acceptor one via the electromag¬ 
netic interaction. The donor and acceptor systems may 
be realized as quantum dots and quantum wells, 1 3 bi¬ 
ological moleculesdefects in semiconductor.® 7 Typi¬ 
cally, the range of the Forster interaction is on the order 
of several nm. 8 

One can try to control the efficiency of the trans¬ 
fer by embedding the donors and acceptors into the 
structured electromagnetic environment. Contradicting 
experiment al rep orts of the enhancecP, modulatecP®, 
independent 1112 and suppressed? 3 transfer are available 
in literature. In the particular case of Ref.[l3]the transfer 
has been studied for the dye molecules embedded in the 
thin polymer film. The speedup of the donor emission de¬ 
cay in the presence of the acceptors has been attributed 
to the energy transfer rate. The transfer rate has been 
suppressed when the film was put on top of the metallic 
mirror. As has been pointed out in Ref. [12j for the pla¬ 
nar mirror this effect cannot be explained by the existing 
well-developed theory of energy transfer between indi¬ 
vidual independent pairs of donors and acceptors! 14 ®^It 
has been suggested in Ref. |13j that the observed trans¬ 
fer suppression can be due to the collective effects, since 
the concentration of the donor dye molecules was quite 
large and the intermolecular distance was on the order 
of several nm. However, any further theoretical expla¬ 
nation was missing in Ref. |l3j Developing a theory of 
collective energy transfer from donors to acceptors in the 
presence of a planar mirror is the main goal of the current 
study. Recently, similar theory has been put forward by 
Pustovit, Urbas and Shahbazyan in Ref.[l8|for molecules 
surrounding the spherical metallic particles. While the 
general approach is similar, the electromagnetic modes 
for the spherical geometry are quite different, so the re¬ 
sults of Ref. m cannot be applied to our system directly. 
Hence, there is a growing demand from both theory and 
experiment to elaborate on the energy transfer in planar 
plasmonic systems with large molecule concentration. 



FIG. 1: Schematic illustration of the studied system. Arrays 
of donor (D) and acceptor (A) molecules above the substrate 
are indicated. 


In the rest of the paper we present the model (Sec. 0 
for the collective transfer and show, that the transfer rate 
can be strongly modified by the mirror (Sec. |III| ). 


II. MODEL 

The system under consideration is schematically de¬ 
picted in Fig. [l] It consists of the metallic substrate 
with the permittivity s su b s bounded from top with the 
semi-infinite dielectric matrix with the permittivity 
Arrays of donor and acceptor molecules form square lat¬ 
tices with the period a, that are parallel to the substrate. 
The donor array is placed at the height h from the sub¬ 
strate, the acceptor array lies in the plane z = h + Sh at 
the vertical distance Sh from the acceptor array. The x 
and y axes are parallel to the interface plane z = 0. 

We use the semiclassical coupled dipole approximation 
with the random source terms, similarly to the approach 
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in Refs. mm 

—PD,j = fD,j + 53 G ( r D,j ~ r D,j')PDJ' , ( 1 ) 

a D 7T. 

^~PA,j = 53 G ( rA ’i ~ r D,j')PD,j' ■ (2) 

a A .. 

3' 

(3) 


rate Y a of the excited acceptor state is assumed enough 
fast to quench both the acceptor-acceptor coupling and 
possible transfer of the energy back to the donors. 

We will first calculate the population of donors by solv¬ 
ing Eq. 0 . This can be accomplished by expanding the 
donor dipole momenta over the Bloch modes, character¬ 
ized by the in-plane Bloch wave vector k: 


Here ppj (paj) are the dipole momenta of the donors 
(acceptors) located at the points rp,j (r A,j ) and having 
the polarizabilities ap {ola), G(r,r f ) is the tensor Green 
function 

VxVxG(r,r')= ^ s(r)G(r, r')+47r f—^ 16(r— r'), 

_ (4) 

accounting for the presence of the metallic mirrorp° and 
fn, j are the random source terms. We assume incoherent 


PD,j = 53 elkrD ’ k PD,k 
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dk x dk y e ikr ^ k p D , k . (8) 
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Substituting Eq. (|8| into Eq. ([!]) we obtain a system of 
independent equations for the Bloch amplitudes Pd k- 
Solving them and substituting the result back to Eq. (jf]) 
we obtain 


stationary pumping of the donor array, in this case the 
realization-averaged values of fpj are equal to zero. The 
pumping to different molecules can be considered inde¬ 
pendent and is characterized by the following correlator: 

Pdj = 53E elfe(rD "' rDJ,) a Dtk f D j> . 

k j' 

where 

(9) 
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OLD,k = 1 n , 

1 — apCp^k 

(10) 

and 


here S is the effective rate of the generation of excitons 
in donors 1 ^, /i, v — x,y,z and the brackets denote the 

CD,k = y ^ e - lkr D,3G(—rpj) 

(11) 

averaging over the realizations of the random sources. 

j/o 


The polarizabilities have the resonant form, 

is the interaction constant for the planar donor array! 22 ^ 


: 


d 2 

a M 


h[uj M —u — i(r M + To,m)] ’ 


M = D,A , (6) 


where dM are the matrix elements of the dipole moment 
for donors and acceptors, ujm are the resonant frequen¬ 
cies and T are the phenomenological nonradiative decay 
rates. The radiative decay rate To is then described by 
the standard expression 21 


r 0 ,M — 


2u\fd 2 M ^/sb 

3 he 3 


(7) 


Our goal is to determine the stationary population of ac¬ 
ceptors induced as a result of the pumping of the donors 
and the energy transfer from donors to acceptors. Im¬ 
portantly, in Eq. 0 we include the electromagnetic cou¬ 
pling between the donor molecules (second term in the 
right hand side). This coupling leads to the formation 
of the collective states, eigenmodes of Eq. 0 . Impor¬ 
tantly, the collective effects are manifested even for inco¬ 
herent pumping. The only necessary conditions are the 
weak inhomogeneous broadening in the donor array and 
long enough lifetime of the donor states. In the linear-in¬ 
pumping regime the pumping determines only the pop¬ 
ulation of the different eigenmodes, but not their spatial 
structure. In order to simplify the model we neglect the 
electromagnetic coupling between acceptors in Eq. §. 
The reason is that the energy is transferred to the ex¬ 
cited states of the acceptors. The nonradiative decay 


Eq. (§ fully accounts for the hybridization between the 
donor array and the electromagnetic modes of the metal¬ 
lic substrate, including the plasmonic modes. This is 
encoded in the polarizabilities of the Bloch modes ap,ki 
which have the form of the bare donor polarizability ap 
renormalized by the interaction with the mirror. The in¬ 
teraction constant Cn,k can be evaluated either by Flo- 
quer summation 22 or by an Ewald-type summatioiPHl, 
here we resort to the latter procedure. 

The average occupation number of the donor molecules 
can be introduced as the averaged squared amplitude of 
the donor dipole momentum divided by the dipole mo¬ 
mentum matrix element: 


1 

MT 2 


Nd — T7~\2 (IP-dWI 2 ) 


\d 




( 12 ) 

Substituting Eq. 0 into Eq. ( [l2| and using the random 
source correlation function Eq. (pi) we obtain 




Nd = S I P> ]Tra nk a\ 


x D,k^D,k 


(13) 


Now we will assume weak coupling between the donor 
array and the mirror. In this approximation the interac¬ 
tion with the mirror leads only to the modification of the 
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donor lifetime. The complex polarizabilities ctD,k then 
have the form 


= 


d D 


(14) 


where 


Fk — 1 + 


h[cdM — u — + ^fcr 0 ,D)] 

ImC k (u D ), (15) 


2 (wd/c) 3 v / £6 


is the tensor analogue of the Purcell factor, describing 
the lifetime modification. The matrix F & can be diago¬ 
nalized, 


F k ,^u = ^ Vfe )At A/fe,A[r fe X ]At/, (J,,V = x,y,z , (16) 

A=1 

where are the Purcell factors characterizing the 

decay of the Bloch modes with different polarizations. 
There exist three possible polarizations for each value of 
k. Note, that the matrix Vk in Eq. (16) is in general case 
neither Hermitian nor symmetric. After the frequency 
integration in Eq. (13) is carried out, the result for the 


donor population assumes the form 

N D = Sj2r k ,x 

fc, A 

where 


Tk, a = 


2(T D + f k ,xV^[V k ^U) 


(17) 


(18) 


is the effective lifetime of the Bloch mode and the sum¬ 
mation over the dummy index /i is assumed. Eq. (17) 


has a clear physical meaning: the total population of the 
donors is a sum over the populations of different eigen- 
modes. Each eigenmode population is given by the prod¬ 
uct of the pumping rate and the lifetime. Since donors 
are pumped in an uncorrelated way, the pumping con¬ 
stant is independent of the Bloch vector k. However, the 
lifetimes can differ for different modes. 

Now we proceed to the calculation of the acceptor pop¬ 
ulation. To this end we substitute Eq. ([9| into Eq. (J2| 
and calculate the correlation function 


N a = 


\ d A\ 


;(\PA(t)\ 2 ) 


ihf I I paAJ)) - (19) 


This yields (cf. with Eq. (13) and see also Ref. [24]) 


Na = S I ^ \a A 




27r 


|2 £ 

k 


Tr Guam ol 


D,k a D,k^DA,k 


G 


t 


where 


GoA,k = y ^e lkrA ’ j G(rAj - r D j 0 


( 21 ) 


is the Green function Fourier component describing the 
donor-acceptor coupling. In the chosen approximation 
the lifetime of the (excited) acceptor states is deter¬ 
mined solely by energy relaxation processes, Fa 
ro,A, r 05 £> 5 T£>, 1/t^\. Hence, we can introduce the effec¬ 
tive energy transfer rate 1 /tet from the following kinetic 
equation for the acceptors: 


Na = Nd 
ta t E t 


( 22 ) 


The left-hand side of this equation describes the (non- 
radiative) decay of acceptors with the lifetime ta = 
1/(2T^), the right-hand-side describes the transfer of en¬ 
ergy from donors with the population N&. Integrating 
over frequency in Eq. (20), using the donor population 
Ne from Eq. ( p~7| ) and the transfer rate definition Eq. (22) 
we obtain 


1 

tet 


2tt 1 


F a 


-\d A \ 2 \d D \ 2 Q 2 


h 7r (uj d - uj a ) 2 + F\ 
where the effective constant of the transfer rate reads 


(23) 


g> = Z = 
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E Tk„ 

k, X 


F- 
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Tk,x\ G DA,k,\\ 


with 


G DA,k,X — en[GDA,k\nvVk,v X 


(24) 


(25) 


Equations (23)—(25) for the effective energy transfer rate 
constitute the central result of this study. Their struc¬ 


ture is quite clear. The first resonant factor in Eq. (23) 


depends on the spectral detuning between donor and ac¬ 
ceptor molecules. The second factor Q 2 describes the 
electromagnetic coupling between the molecules arrays. 
If the collective coupling between the donors is neglected, 
the factor reduces to 

So 2 = £ Tr [G\r A ,j - r DJ >)G(r A ,j - r DJ >)], (26) 


i.e. it is described by a sum of individual transfer events 
from different donors. In this case the transfer rate as¬ 
sumes the same form as in Refs. EHM In the general 
case, however, the transfer constant Q 2 is given by a sum 
of the contributions from Bloch eigenmodes of the donor 
ensemble. Since the donors are pumped incoherently and 
uncorrelated, the contributions are independent. Each 
term in the sum in Eq. (23) is proportional to the life¬ 


( 20 ) 


time of the corresponding mode r k,x> that determines the 
mode population. 


III. RESULTS AND DISCUSSION 

In this section we present detailed analysis of the ef¬ 
fect of the metallic mirror on the transfer rate Eq. (24). 
Figure [2] shows the dependence of the transfer rate on the 
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FIG. 2: (Color online) Energy transfer rate as function of 
the height of the donors above the mirror h for (a) acceptors 
located above donors and (b) acceptors located in the same 
horizontal plane as donors. Thick solid/red and dashed/blue 
curves correspond to the transfer rates calculated according to 
Eq. (24) with the lattice constants a — 7 nm and a — 15 nm, 
respectively. The curves are normalized to the transfer rates 
between single donor and single acceptor at the corresponding 
distance without the mirror. Thin curves show the transfer 
rates between single donor and single acceptor, i.e. calculated 
from Eq. (26) neglecting the collective effects. The calculation 
parameters are as follows: A = 2ttc/ujd = 500 nm, £b = 2, 
Td = 0.1To,d, £ S ubs = —9.77 + 0.31i (corresponds to Ag in 


Ref. l25l) , vaj — vdj = A hz, Ah — 3.5 nm (a) and vaj — 
ro,j — a/4(x + y) (b). The insets schematically illustrate the 
geometry of the problem. 


height of the donor and acceptor molecules above the mir¬ 
ror. The lattice constant of the donor and acceptor arrays 
and their relative positions with respect to each other re¬ 
main fixed. Two panels of Fig. [ 2 ] correspond to different 
geometries: (a) out-of-plane transfer, when acceptors are 
positioned above the donors at the height 5h and (b) in¬ 
plane transfer, when acceptors are positioned in the same 
horizontal plane as the donors so that each acceptor lies 
in the center of the square unit cell of the donor lattice. 
Thick lines are calculated according to Eq. (25) and fully 
account for the hybridization between the donors while 
thin lines are calculated in the independent donor ap¬ 
proximation Eq. (26) by retaining only one term for one 
donor, that is closest to the acceptor. Importantly, the 
transfer rates in the collective model and the indepen¬ 
dent donor model are different from each other even in 
the absence of the mirror (or for h a). The reason is 
that the local field of the donor array at the acceptor is 
different from that of the single donor even in vacuum. In 
order to compare the two models visually, the curves for 
the single-donor model have been scaled by the factors 
0.3 and 1.4 in Fig. if and Fig. if, respectively. The main 
result of this work, shown in Fig. [2j is as follows: the col¬ 
lective transfer rate is strongly sensitive to the distance 
from the mirror (thick lines) while the independent-donor 
one is not. When the acceptors are positioned above the 
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FIG. 3: Purcell factors for the Bloch modes of the donor ar¬ 
ray with different polarizations as functions of the Bloch wave 
vector along the T-X direction. Thick/black and thin/red 
curves correspond to h— 10 nm and h — 20 nm, respectively. 
Calculated for a — 7 nm and the same other parameters as in 
Fig. [2] The arrow indicates the value of k equal to the light 
cone edge uJDy/sb/c. 


donors, the transfer is strongly suppressed by the mir¬ 
ror, while for acceptors located between the donors the 
transfer is enhanced. This effect cannot be described in 
the independent transfer model (thin curves), which pre¬ 
dicts almost vanishing dependence of the transfer rate on 
the distance up to the distance h ~ 3 nm, in agreement 
with Refs. 1121131 The transfer modification is most pro¬ 
nounced when the height of the donor array above the 
mirror is on the order of the array period a (cf. solid 
and dotted curves for a = 7 nm and a = 15 nm). As the 
period increases, the effect is suppressed. Thus, collec¬ 
tive transfer model predicts a strong modification of the 
transfer rate near the metallic mirror, in stark contrast 
to the independent transfer model. Now we proceed to 
the explanation of these findings. 

Our analysis indicates, that the Green functions, enter¬ 
ing Eq. (25), are not very sensitive to the distance from 
the mirror. Hence, the strong transfer modification, ob¬ 
served in Fig. [2] is solely due to the dependence of the 
lifetimes of the Bloch eigenmodes Tk,\ on the distance. 
The Purcell factors fk,\, characterizing the lifetime mod¬ 
ification by the mirror, are presented in Fig. [3] Thick 
and thin curves were calculated as functions of the Bloch 
wave vector for two distances, h = 10 nm and h = 20 nm, 
respectively. For each height we show three curves, corre¬ 
sponding to different eigenmode polarizations. The cal¬ 
culation reveals strong dependence of the Purcell factors 
on the wave vector, distance and polarization. For given 
distance one can distinguish between several regions in 
Fig. [3j The first region is inside the light cone {k <u/c). 
In this case the Purcell factor is much larger than unity 
even at the large distances from the mirror (or without 
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the mirror). This can be understood as a semiclassi- 
cal counterpart of the Dicke superradiancepS when the 
donors, that are separated by the distances smaller than 
the light wave length, emit collectively. The radiative 
rate is larger by the factor on the order of l/(ua/c ) 2 
than that of the individual donor in vacuum. 27 ' Second, 
one can see the peak in Fig. [3] at /c ~ 0.047r/a, that de¬ 
scribes strong enhancement of the donor decay due to the 
plasmons. The peak is present only for one of the curves, 
corresponding to the transverse magnetic (TM, or p) po¬ 
larization. Third, there exists a region of larger Bloch 
wave vectors k > 0.l7r/a. This range corresponds to the 
evanescent modes, with the electric field exponentially 
decaying from the donor plane and lying outside the light 
cone. In the absence of the mirror these waves would be 
completely evanescent and they would have zero radiative 
decay rate, /& = 0. Due to the presence of the mirror, 
the modes acquire finite lifetime due to the absorption 
in the substrate, /*, oc Im5 su b s e _2fc/l P- The weakest val¬ 
ues of the Purcell factor in this third region correspond 
to the donor modes with the polarization close to trans¬ 
verse electric (TE, or 5 ). Although the Purcell factor 
for the evanescent modes is much less than unity, it is 
strongly enhanced when the array approaches the mirror 
and the coupling to the mirror increases (cf. black and 
red curves). We will now explain how this very effect 
leads to the modification of the transfer rate Eq. (24) by 
the mirror. 

Figure [4] shows the dependence of the Green func¬ 
tion components Eq. ( [25] ) and the lifetimes of donor 
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Bloch wave vector 
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FIG. 4: Green function Eq. (25) characterizing the transfer 
strength (a) and lifetimes of the Bloch modes of the donor 
array (b) as functions of the Bloch wave vector in the T— 
M direction [inset in panel (b)]. Black/solid and red/dotted 
curves correspond to the out-of-plane energy transfer and in¬ 
plane transfer, as indicated by the insets in panel (a). For each 
geometry only the mode with one polarization, providing the 
largest contribution to the energy transfer, is shown. Other 
calculation parameters as the same as in Fig. [2] 


modes on the Bloch wave vector for out-of-plane (black 
solid curves) and in-plane transfer (red dotted curves). 
For each geometry we have plotted the Green function 
component only for one polarization, that is largest by 
the absolute value and provides the dominant contribu¬ 
tion to the transfer rate. The Green function values in 
panel (a) have been additionally multiplied by the fac¬ 
tor &, accounting for the growth of the density of states 
d 2 k = kdkdifk with the absolute value of the wave vec¬ 
tor. The transfer rate Eq. ( [24] ) is proportional to the 
integral of the product of the values of k\GDA,k,x\ 2 i n 
Fig. [4](a) and t^a in Fig. |4^b) over the Bloch wave vec¬ 
tor k. Importantly, due to the difference in the Green 
function values \GDA,k,x\ 2 , the out-of-plane transfer is 
mediated by the p-polarized Bloch modes, while the in¬ 
plane transfer is controlled by the 5-like modes. These 
modes have quite different lifetimes, that are plotted in 
Fig# for the height h — 20 nm above the mirror (see 
also Fig. [2|. This explains the qualitative difference be¬ 
tween the dependence of the in-plane transfer and the 
out-of-plane transfer on the distance from the mirror. In 
order to elucidate the effect we will now analyze these two 
cases in more detail, (i) Out-of-plane transfer. When the 
array approaches the mirror, the lifetime of the p-type 
modes, determining the transfer, becomes shorter. The 
product W = Y.k,\ T k,x\ G DA,k,\\ 2 , controlled by the p- 
type modes, is decreased. On the other hand, the total 
lifetime T = A t^a of the donor array is due to the 
long-living modes with 5-type polarization, rather than 
the p modes. The lifetime of 5-type modes is equal to 
1/2Td for k > 0.l7r/a and weakly sensitive to the mirror 
(red/dotted curve in Fig. [ 4 ) 3 ). As a result, W decreases 
faster than the T when the mirror is brought closer, and 
the transfer rate Eq. ( [24] ) oc W/T is suppressed, in agree¬ 
ment with Fig. [2ja) . (n) In-plane transfer. In this case 
the function W is determined by the long-living 5-type 
eigenmode. Hence, the value of W is practically insensi¬ 
tive to the mirror because the lifetime is modified only for 
very small k < 0.l7r/a, where the product of the squared 
Green function and the density of states is small. Hence, 
W decreases slower than T, and the transfer rate W/T is 
enhanced at smaller distances from the mirror, in agree¬ 
ment with Fig. [2jb) . 


Finally, in Fig. [5] we analyze the dependence of the 
transfer rates on the nonradiative decay of the donors 
r^. Clearly, when the nonradiative component of the 
decay rate increases, the decay rate of any given Bloch 
mode becomes less sensitive to the mirror presence. As 
a result, the total energy transfer rate returns to its bulk 
value, unaffected by the mirror. This is in agreement with 
the result of calculation in Fig. [5] both for out-of-plane 
(red/dotted curve) and in-plane (black/solid) geometries 
the energy transfer tends to its bulk value for large Td- 
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FIG. 5: (Color online) Energy transfer rate as function of 
the nonradiative decay of the single donor T. Black/solid 
and red/dotted curves correspond to the out-of-plane energy 
transfer and in-plane transfer, as indicated by the insets. 
Other calculation parameters as the same as in Fig. [2] 

IV. SUMMARY 

To summarize, we have developed a theory of energy 
transfer between the planar arrays of donor and acceptor 
molecules, lying on top of a metallic mirror in the regime 
of weak, incoherent and spatially uncorrelated pumping 
of the donors. We have used a coupled dipole model with 
random source terms that has allowed us to obtain ana¬ 
lytical results for the transfer rate. It has been demon¬ 
strated, that the transfer rate between the molecule ar¬ 


rays can be strongly modified by the mirror, in stark con¬ 
trast to the transfer between individual molecules. The 
enhanced sensitivity of the transfer to the presence of the 
mirror is due to the modification of the lifetimes of the 
collective Bloch modes of the donor array. Depending 
on the relative position between the donor and acceptor 
molecules, the transfer can be both suppressed (accep¬ 
tors between the donors) or enhanced (acceptors above 
the donors). While further study is needed to account 
for the possible effect of the strong coupling between the 
donors and the plasmonsj^ for inevitable effects of dis¬ 
order in actual samples, and for the bulk film geometry, 
these our results might already provide the key to ex¬ 
plain the suppression of the energy transfer in the plas- 
monic environment, recently observed experimentally in 
Ref. H3J 
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